#install.packages("measurements")
library(ggmap)
library(measurements)
pcs <- read_sf("alldata/Scotland_pcs_2011")
#Print Post codes lists
arransubsect <- filter(pcs,substr(label,1,4)=="KA27")
arransubsect %>%
  ggplot() +
  geom_sf() +
  theme(axis.text.x=element_text(angle=45, hjust = 1))

head(arransubsect)
Simple feature collection with 1 feature and 2 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 186130 ymin: 618830.5 xmax: 206925.5 ymax: 652393
epsg (SRID):    NA
proj4string:    +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs
    name  label                       geometry
1 KA27 8 KA27 8 MULTIPOLYGON(((188338 63101...
ggplot() +
  geom_point(data= arrancoordinates, mapping = aes(x = longitude, y = latitude))

ggplot() +
  geom_point(data= arrancoord2, mapping = aes(x = longitude, y = latitude))

head(arrancoordinates)
arrancoord2<- arrancoordinates
arrancoord3<- arrancoordinates
library("sp")
?char2dms
# change the space symbol to a degree
arrancoord2$latitude = gsub(' ', '°', arrancoord2$latitude)
arrancoord2$longitude = gsub(' ', '°', arrancoord2$longitude)
# convert from decimal degrees to decimal minutes
arrancoord2$latitude = measurements::conv_unit(arrancoord2$latitude, from = 'dec_deg', to = 'deg_dec_min')
arrancoord2$longitude = measurements::conv_unit(arrancoord2$longitude, from = 'dec_deg', to = 'deg_dec_min')
head(arrancoord2)
coordinates <- select(arrancoordinates, latitude, longitude)
coordinates
arrancoordinates %>% 
    select(latitude, longitude) %>% 
    head
# change the degree symbol to a space
#arrancoordinates$latitude = gsub(' ', '°', arrancoordinates$latitude)
#arrancoordinates$longitude = gsub(' ', '°', arrancoordinates$longitude)
# convert from decimal minutes to decimal degrees
#arrancoordinates$latitude = #measurements::conv_unit(arrancoordinates$latitude, from = 'dec_deg', to = 'deg_dec_min')
#arrancoordinates$longitude = #measurements::conv_unit(arrancoordinates$longitude, from = 'dec_deg', to = 'deg_dec_min')
msleep %>% 
    mutate(rem_proportion = sleep_rem / sleep_total) %>%
    head
msleep %>% 
    mutate(rem_proportion = sleep_rem / sleep_total, 
           bodywt_grams = bodywt * 1000) %>%
    head
require(maps)
d = data.frame(map(database="italy", plot=F)[c("x", "y")])
ggplot() +
coord_map() +
geom_polygon(data=d, mapping=aes(x=x, y=y))

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CiNpbnN0YWxsLnBhY2thZ2VzKCJtZWFzdXJlbWVudHMiKQpgYGAKCgpgYGB7cn0KbGlicmFyeShnZ21hcCkKbGlicmFyeShtZWFzdXJlbWVudHMpCmBgYAoKCmBgYHtyfQpwY3MgPC0gcmVhZF9zZigiYWxsZGF0YS9TY290bGFuZF9wY3NfMjAxMSIpCgojUHJpbnQgUG9zdCBjb2RlcyBsaXN0cwphcnJhbnN1YnNlY3QgPC0gZmlsdGVyKHBjcyxzdWJzdHIobGFiZWwsMSw0KT09IktBMjciKQphcnJhbnN1YnNlY3QgJT4lCiAgZ2dwbG90KCkgKwogIGdlb21fc2YoKSArCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdCA9IDEpKQpgYGAKCmBgYHtyfQpoZWFkKGFycmFuc3Vic2VjdCkKYGBgCgoKYGBge3J9CmdncGxvdCgpICsKICBnZW9tX3BvaW50KGRhdGE9IGFycmFuY29vcmRpbmF0ZXMsIG1hcHBpbmcgPSBhZXMoeCA9IGxvbmdpdHVkZSwgeSA9IGxhdGl0dWRlKSkKYGBgCmBgYHtyfQpnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhPSBhcnJhbmNvb3JkMiwgbWFwcGluZyA9IGFlcyh4ID0gbG9uZ2l0dWRlLCB5ID0gbGF0aXR1ZGUpKQpgYGAKCmBgYHtyfQpoZWFkKGFycmFuY29vcmRpbmF0ZXMpCmBgYAoKYGBge3J9CmFycmFuY29vcmQyPC0gYXJyYW5jb29yZGluYXRlcwpgYGAKCmBgYHtyfQphcnJhbmNvb3JkMzwtIGFycmFuY29vcmRpbmF0ZXMKYGBgCgpgYGB7cn0KbGlicmFyeSgic3AiKQo/Y2hhcjJkbXMKYGBgCgpgYGB7cn0KCiMgY2hhbmdlIHRoZSBzcGFjZSBzeW1ib2wgdG8gYSBkZWdyZWUKYXJyYW5jb29yZDIkbGF0aXR1ZGUgPSBnc3ViKCcgJywgJ8KwJywgYXJyYW5jb29yZDIkbGF0aXR1ZGUpCmFycmFuY29vcmQyJGxvbmdpdHVkZSA9IGdzdWIoJyAnLCAnwrAnLCBhcnJhbmNvb3JkMiRsb25naXR1ZGUpCgojIGNvbnZlcnQgZnJvbSBkZWNpbWFsIGRlZ3JlZXMgdG8gZGVjaW1hbCBtaW51dGVzCmFycmFuY29vcmQyJGxhdGl0dWRlID0gbWVhc3VyZW1lbnRzOjpjb252X3VuaXQoYXJyYW5jb29yZDIkbGF0aXR1ZGUsIGZyb20gPSAnZGVjX2RlZycsIHRvID0gJ2RlZ19kZWNfbWluJykKYXJyYW5jb29yZDIkbG9uZ2l0dWRlID0gbWVhc3VyZW1lbnRzOjpjb252X3VuaXQoYXJyYW5jb29yZDIkbG9uZ2l0dWRlLCBmcm9tID0gJ2RlY19kZWcnLCB0byA9ICdkZWdfZGVjX21pbicpCmBgYAoKYGBge3J9CmhlYWQoYXJyYW5jb29yZDIpCmBgYAoKCmBgYHtyfQpjb29yZGluYXRlcyA8LSBzZWxlY3QoYXJyYW5jb29yZGluYXRlcywgbGF0aXR1ZGUsIGxvbmdpdHVkZSkKYGBgCgpgYGB7cn0KY29vcmRpbmF0ZXMKYGBgCgpgYGB7cn0KYXJyYW5jb29yZGluYXRlcyAlPiUgCiAgICBzZWxlY3QobGF0aXR1ZGUsIGxvbmdpdHVkZSkgJT4lIAogICAgaGVhZApgYGAKCmBgYHtyfQojIGNoYW5nZSB0aGUgZGVncmVlIHN5bWJvbCB0byBhIHNwYWNlCiNhcnJhbmNvb3JkaW5hdGVzJGxhdGl0dWRlID0gZ3N1YignICcsICfCsCcsIGFycmFuY29vcmRpbmF0ZXMkbGF0aXR1ZGUpCiNhcnJhbmNvb3JkaW5hdGVzJGxvbmdpdHVkZSA9IGdzdWIoJyAnLCAnwrAnLCBhcnJhbmNvb3JkaW5hdGVzJGxvbmdpdHVkZSkKCiMgY29udmVydCBmcm9tIGRlY2ltYWwgbWludXRlcyB0byBkZWNpbWFsIGRlZ3JlZXMKI2FycmFuY29vcmRpbmF0ZXMkbGF0aXR1ZGUgPSAjbWVhc3VyZW1lbnRzOjpjb252X3VuaXQoYXJyYW5jb29yZGluYXRlcyRsYXRpdHVkZSwgZnJvbSA9ICdkZWNfZGVnJywgdG8gPSAnZGVnX2RlY19taW4nKQojYXJyYW5jb29yZGluYXRlcyRsb25naXR1ZGUgPSAjbWVhc3VyZW1lbnRzOjpjb252X3VuaXQoYXJyYW5jb29yZGluYXRlcyRsb25naXR1ZGUsIGZyb20gPSAnZGVjX2RlZycsIHRvID0gJ2RlZ19kZWNfbWluJykKYGBgCgpgYGB7cn0KCmBgYAoKCmBgYHtyfQptc2xlZXAgJT4lIAogICAgbXV0YXRlKHJlbV9wcm9wb3J0aW9uID0gc2xlZXBfcmVtIC8gc2xlZXBfdG90YWwpICU+JQogICAgaGVhZApgYGAKCgpgYGB7cn0KbXNsZWVwICU+JSAKICAgIG11dGF0ZShyZW1fcHJvcG9ydGlvbiA9IHNsZWVwX3JlbSAvIHNsZWVwX3RvdGFsLCAKICAgICAgICAgICBib2R5d3RfZ3JhbXMgPSBib2R5d3QgKiAxMDAwKSAlPiUKICAgIGhlYWQKYGBgCgoKYGBge3J9CnJlcXVpcmUobWFwcykKZCA9IGRhdGEuZnJhbWUobWFwKGRhdGFiYXNlPSJpdGFseSIsIHBsb3Q9RilbYygieCIsICJ5IildKQpgYGAKCmBgYHtyfQpnZ3Bsb3QoKSArCmNvb3JkX21hcCgpICsKZ2VvbV9wb2x5Z29uKGRhdGE9ZCwgbWFwcGluZz1hZXMoeD14LCB5PXkpKQpgYGAKCg==